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1.  Introduction 


The  Advanced  Munitions  Concepts  Branch  of  the  U.S.  Army  Research  Laboratory’s  (ARL’s) 
Weapons  and  Materials  Research  Directorate  has  for  many  years  designed,  built,  and  employed 
body-fixed  sensor  and  telemetry  systems  to  measure  flight  body  kinematics,  primarily  in  support  of 
military  ordnance  testing.  Because  requirements  imposed  by  military  applications  such  as  high-g 
environment,  extreme  projectile  dynamics,  small  size,  low  cost,  low  power  consumption,  etc., 
exclude  many  traditional  inertial  sensor  systems,  ARL  is  continually  exploring  emerging  tech¬ 
nologies  and  developing  alternate  techniques  for  their  utility  in  obtaining  desired  measurements. 

With  the  result  from  vector  differentiation  relating  the  time  derivatives  of  a  vector  represented  in 
two  coordinate  systems  in  relative  motion,  a  technique  for  obtaining  pitch,  yaw,  and  roll  rates  from 
a  single,  low-cost,  commercial  off-the-shelf  magnetometer  was  developed.  In  a  recent  ARL  flight 
experiment  intended  to  characterize  the  angular  motion  of  National  Aeronautics  and  Space  Admin¬ 
istration  (NASA)  re-entry  capsules,  no  data  were  available  from  the  strap-down  angular  rate 
sensors  during  critical  portions  of  several  flights  because  the  launch  environment  and  projectile 
dynamics  exceeded  sensor  capabilities.  These  rates  were  successfully  derived  from  magnetometer 
data,  and  complete  Euler  angular  histories  of  the  test  trajectories  were  obtained  with  an  ARL- 
developed  vector-matching  algorithm. 

In  this  report,  the  magnetometer-based  methodology  is  presented,  the  flight  experiment  and  sub¬ 
sequent  analyses  are  described,  criteria  for  use  of  this  methodology  are  given,  and  potential  uses  of 
this  technique  in  inertial  measurement  unit  (IMU)/inertial  navigation  system  (INS)  applications  are 
discussed. 


2.  Inertial  Navigation  of  Projectiles 


The  equations  of  motion  describing  free-flight  dynamics  of  rigid  projectiles  have  six  degrees  of 
freedom,  three  translational  velocity  components  and  three  rotational  velocity  components.  A 
traditional  strap-down  IMU  consists  of  an  orthogonal  triad  of  linear  accelerometers  and  an  or¬ 
thogonal  triad  of  angular  rate  sensors  oriented  along  a  projectile’s  principal  axes.  Given  initial 
launch  position,  orientation,  and  velocity  of  a  projectile,  the  rate  sensors’  output  is  integrated  to 
update  the  projectile  orientation,  and  the  accelerometers  are  integrated  once  to  update  the  projectile 
velocity  and  twice  to  update  the  projectile  position.  The  solution  of  the  inertial  navigation  problem 
is  conceptually  simple,  but  it  is  often  difficult  to  realize  an  IMU  capable  of  measuring  the  required 
six  body  states  under  constraints  imposed  by  a  particular  application.  Although  this  problem  is 
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equally  of  concern  to  accelerometers  and  rate  sensors,  only  the  orientation  estimation  (i.e.,  rate 
sensor)  problem  is  addressed  herein. 

Formulation  of  the  inertial  navigation  problem  for  gun-  and  tube-launched  projectiles  requires  the 
use  of  multiple  coordinate  systems  (Flarkins,  2003,  2007).  Trajectory  time  histories  are  best 
described  in  an  earth-fixed  coordinate  system  with  its  origin  at  the  launcher.  Of  necessity,  strap- 
down  sensor  measurements  are  made  in  a  flight-body-fixed  coordinate  system,  and  target  locations 
are  most  naturally  described  in  another  earth-fixed  system. 

The  first  coordinate  system  is  right-handed  Cartesian  (I,J,K)  with  its  origin  at  the  launch  site.  This 
will  be  referred  to  as  the  “earth-fixed”  system  and  the  axes  are  defined  by 

•  The  /  and  J  axes,  which  define  a  plane  tangent  to  the  earth’s  surface  at  the  origin; 

•  The  K  axis,  which  is  perpendicular  to  the  earth’s  surface  with  positive  downward,  i.e.,  in 
the  direction  of  gravity; 

•  The  /  axis,  which  is  chosen  so  that  the  centerline  of  the  launcher  is  in  the  /-if  plane. 

Down-range  travel  is  then  measured  along  the  /  axis,  deflection  along  the  J  axis  (positive  to  the 
right  when  one  is  looking  down  range),  and  altitude  along  the  K  axis  (positive  downwards)  (see 
figure  1). 


The  second  system  is  convenient  for  aeroballistic  computations  of  rigid  projectiles’  flights  and  for 
describing  the  locations  and  orientations  of  such  projectiles’  components.  This  system  is  right- 
handed  Cartesian  (i,j,k)  with  its  origin  at  the  center  of  gravity  (c.g.)  of  the  flight  body.  For  rotating 
flight  bodies,  the  projectile-fixed  coordinate  system  usually  has  its  i  axis  lying  along  the  projectile 
axis  of  symmetry,  i.e.,  the  spin  axis  (with  positive  in  the  direction  of  travel  at  launch).  The  j  and  k 
axes  are  then  oriented  so  as  to  complete  the  right-handed  orthogonal  system  (figure  1).  Spin  (p ), 
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pitch  ( q ),  and  yaw  (r)  rates  are  measured  about  these  axes.  This  will  be  referred  to  as  the  “body- 
fixed”  system. 

The  third  coordinate  system  ( X,  Y,Z)  is  commonly  employed  to  specify  locations  on  or  near  the 
earth’s  surface,  i.e.,  north,  east,  and  down.  This  will  be  referred  to  as  the  “navigation”  system 
where  north  =  X,  east  =  Y,  and  down  =  Z. 


The  earth-fixed  and  body-fixed  coordinate  systems  are  related  through  an  Euler  rotation  sequence, 
beginning  with  a  rotation  of  the  earth-fixed  frame  about  the  K-axis  through  the  yaw  angle  y/  .  The 

system  is  then  rotated  about  the  new  J’-axis  through  the  pitch  angle  6  .  Finally,  the  system  is 
rotated  about  the  new  i-axis  through  the  roll  angle  (f) .  The  two  systems  are  related  by  the  direction 

cosine  transformation  matrix  (DCM),  TEb  ,  with  the  subscript  denoting  earth  fixed  to  body  fixed. 

This  transformation  matrix  is 


TEb  ~ 


c¥ce 

cy/s6S(j)  ~'sy/C(j> 
ycy/sOc<j>  sy/s (j> 


sy/c6 

sy/s6s(j)  ~>rCy/c(j) 
s\f/sOc(j>  —  Cy/S(j> 


~se 

C6S(j)  > 

C6C(, j)  j 


(1) 


where  c,  is  cos(») ,  and  s,  is  s i n (• ) .  Figure  2  shows  both  coordinate  systems  and  the  Euler  angle 
relations  between  them. 


Figure  2.  Earth-fixed  and  body-fixed  systems 
and  the  Euler  angle  rotation. 


The  body-fixed  components  of  projectile  angular  motion  and  the  Euler  angle  derivatives  are 
related  by 

<j>  =  p  +  [q  sin(^)  +  r  cos(^)  ]tan(#) 

6  =  q  cos(^)  -  r  sin(^)  (2) 

\f/  =  (q  sin(^)  +  r  cos(tf>))/cos(0 ) 
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3.  Vector  Magnetometer 


Among  the  many  varieties  of  magnetic  sensors,  “vector”  magnetometers  are  devices  whose  output 
is  proportional  to  the  magnetic  field  strength  along  the  sensor’s  axis(es).  If  a  tri-axial  vector  mag¬ 
netometer  is  installed  so  that  the  sensor  axes  are  parallel  to  the  axes  of  the  body-fixed  system,  the 
projections  of  the  earth’s  magnetic  field  onto  each  of  the  sensor  axes  can  be  obtained  by  equation  1 
If  Me={Mi,Mj,MK)  is  the  magnetic  field  vector  in  the  earth-fixed  system,  then  the 
components  along  the  sensor  axes  are  given  by 

Mb  =  TEb  Me  (3) 


or 


Mi  =cy/c0MI  +sy/c0MJ  ~SQMK 


ipy,ses4  -SyC^l  +(?y,sesf  +Cy/C(j)Mj  +CQS(j)M K 
ipy/s0C(f)  sy/s(j)^^ I  ~>r^y/S0C<f>  c9c(j)^ K 


(4) 


In  any  real  magnetic  sensor,  detennination  of  axes’  orientations  and  calibration  coefficients  can  be 
a  complex  process,  but  for  the  present  purpose,  it  is  assumed  that  this  has  been  successfully  accom¬ 
plished  and  Mb  is  being  accurately  measured  at  a  known  sampling  rate. 


4.  Obtaining  Angular  Rates  From  Vector  Magnetometers 


Consider  two  coordinate  systems  with  the  same  origin  and  in  relative  motion,  e.g.,  the  earth- fixed 
and  body-fixed  systems  just  described  at  the  time  of  projectile  launch.  From  vector  differentiation, 
the  time  derivative  of  any  vector  in  the  earth- fixed  system  (vE  =  8  vE  /  8 1)  and  its  time  derivative 
in  the  body-fixed  system  (vb  =  Svb  / 8 1 )  are  related  by 

vE=Vb+a>bxvb  (5) 

where  cob  =  ( p,q,r )  .  Applied  to  the  geomagnetic  field  vector,  equation  5  becomes 

ME=Mb+cdbxMb  (6) 

Realizing  that  equation  6  is  unaffected  by  a  translation  of  the  earth- fixed  system’s  origin  to  the 
projectile  e.g.  at  each  sampling  time  and  that  M E  is  unchanging  in  the  earth-fixed  system  and 
expanding  in  component  form,  we  have  the  relations 

Mj  =  ~qMk  +rM j,  M j  =  pMk  -rMh  and  Mb  =  -pMj  +qMt.  (7) 
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Because  these  equations  are  not  linearly  independent,  they  can  not  be  solved  directly  for  the 
angular  rates.  However,  for  most  rolling  projectiles  where  \p\  »  \q\  and  |r| ,  a  good  estimate  of p 

is  readily  obtainable  from  the  magnetometer  data,  as  described  in  the  next  section.  With  a  spin 
estimate,  this  system  can  be  solved  to  yield  estimates  of  the  body-fixed  pitch  and  yaw  rates. 
Therefore, 

q  =  (Mk+  pMj  )/ Mt  and  r  =  (-  Mj  +  pMk )/ M ,  .  (8) 


5.  Obtaining  Spin  Rate  From  Magnetometers 

Consider  an  earth-fixed,  right-handed  Cartesian  coordinate  system  where  the  z-axis  is  along  the 
geomagnetic  field.  In  this  new  system,  denoted  by  the  subscript  m,  the  geomagnetic  field  vector  is 
Mm  =  (0,0,  M  p-  )  .As  seen  in  section  2,  there  is  a  new  set  of  Euler  angles  that  defines  a  trans¬ 
formation  matrix  from  this  magnetic  coordinate  system  into  the  body-fixed  system  so  that 

Ya\  (  \ 

Cy,  CQ  Sy,  Cg  -Sg  0  -Sg 

T  m  m  i  m  m  m  m 

cw  s0  s<f  ~sw  C(j>  sw  s0  s <b  +  cw  c(j>  c0  s 6  0  —  Cg  Sa  (9) 

t  m  m  t  m  T  m  T  m  r  m  m  t  m  r  m  r  m  m  t  m  m  t  m 

cw  s0  cih  sy/  sd>  sw  s0  cd>  ~~cy/  sd>  c6  cd>  c0  c(b  , 

t  m  m  t  m  r  m  r  m  T  m  m  t  m  t  m  t  m  m  t  m  /  x  ^  \  m  r  m  / 


This  gives  a  definition  of  the  magnetometer  measurements  in  terms  of  the  magnetic  Euler  angles. 
The  magnetic  pitch  angle, 

0m=sin-1(-Mi/\ME\),  (10) 

is  the  complement  of  the  angle  between  the  projectile’s  spin  axis  ( /  ),  and  the  magnetic  field,  M E  . 
The  magnetic  roll  angle,  (j)m  ,  is  computed  by 

<f>m  =tan -l{Mj/Mk)  (11) 

Analogous  to  equation  2,  the  body-fixed  rates  and  the  derivatives  of  the  magnetic  Euler  angles  are 
related  by 

(f>m  =P  +  [q sin(^m )  +  r  cos(^H )  ] tan(6»w ) 

dm  =  q  cos  ((j)m  )  -  r  sin  (<f>m  )  (12) 

Vm  =  (q  sin  (K  )  +  r  cos  (</>m  ))/cos(6»w  ) 
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Estimates  of  <j)m  can  be  obtained  in  several  ways.  The  simplest  method  is  to  make  roll  period 

estimates  from  successive  zero  crossings  or  signal  extrema  on  the  Mj  or  Ma  signals.  This  process 
yields  average  roll  rates  over  the  respective  periods.  More  continuous,  higher  order  estimates  are 
obtained  by  the  computation  of  equation  1 1  at  each  sampling  time  and  the  differentiation  of  the 
results.  Alternatively,  (j))n  is  computed  directly  from 


S tan-1  (Mj/Mk)  _S(Mj/Mk) 

r  i  ' 

= 

(  M jMk  -MkM  j  ^ 

ST  ~  ST 

Kl  +  (Mj/Mk)2  ^ 

{  Mj2+Mk2  J 

with  the  advantage  of  avoiding  potential  singularities  in  equation  1 1  when  =  0  .  The  spin  rate 
( p )  can  then  be  estimated  by  low-pass  fdtering  of  the  (f>m  estimates  (Wilson,  2004). 


6.  Measuring  Angular  Rates  of  a  NASA  Crew  Exploration  Vehicle  (CEV) 
Model 


NASA  needs  to  characterize  the  aerodynamics  of  the  CEV  that  will  be  a  part  of  future  Mars 
missions.  Some  previous  measurements  had  been  made  in  spark  ranges  with  scale  models  of  the 
CEV,  but  this  methodology  cannot  be  employed  to  characterize  all  conditions  of  interest  because  of 
velocity  and  stability  limitations  imposed  by  safety  considerations  in  an  indoor  range.  Further, 
only  limited  amounts  of  data  are  collected  for  each  shot  in  a  spark  range,  so  testing  costs  quickly 
mount  with  the  number  of  shots  required.  With  the  dual  hope  of  expanding  the  set  of  potentially 
measurable  flight  dynamics  and  reducing  testing  costs,  it  was  decided  that  gun  launching  of  Mars 
CEV  models  equipped  with  a  sensor  and  telemetry  system  (figure  3  c)  at  an  outdoor  range  would  be 
explored  as  a  practicable  way  to  acquire  the  desired  data  at  reentry  velocities  (Brown  et  al.,  2006). 
Before  proceeding  to  the  Orion  CEV  tests,  we  evaluated  the  proposed  methodology  using  an 
Apollo  capsule  model  with  known  aerodynamics  (figure  3  a  and  b).  The  sensor  system  consisted  of 
six  angular  rate  sensors  and  a  three-axis  magnetometer.  Along  each  of  the  principal  axes  there 
were  an  angular  rate  sensor  with  a  dynamic  range  of  ±  1000  deg/s,  an  angular  rate  sensor  with  a 
dynamic  range  of  ±  2000  deg/s,  and  a  vector  magnetometer. 

Figure  4  gives  the  body-fixed  pitch  axis  rate  sensor  data  for  the  first  0.5  second  of  one  of  the  Apollo 
model  flight  tests.  Two  “problems”  with  these  data  are  readily  apparent.  First,  the  pitch  (and  yaw) 
angular  rates  exceeded  the  dynamic  range  of  the  rate  sensors  and  clipping  resulted.  Second,  after 
gun  launch,  the  rate  sensors  required  time  to  “settle”.  This  is  obvious  in  the  1000-deg/s  sensor  data 
but  was  later  discovered  to  be  equally  true  of  the  2000-deg/s  sensor  data.  Because  of  these  issues, 
the  magnetometer-based  method  was  used  to  estimate  the  angular  rates  and  the  Apollo  model’s 
attitude  history. 
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Figure  3.  NASA  re-entry  vehicle  models  with  sensor  and  telemetry  system. 


Figure  4.  Apollo  model  flight  data. 

Initial  spin  rate  was  estimated  to  be  approximately  2  Hz  from  a  period  measurement  of  radial 
magnetometer  output.  With  this  value  for  p,  equation  8  was  evaluated  to  obtain  estimates  of  the 

body-fixed  pitch  (y  )and  yaw(f)  rates.  We  estimated  M j  and  Mp  by  differencing  the  successive 

magnetometer  measurements.  The  resulting  pitch  rate  estimate  is  seen  in  figure  5a  superimposed 
on  the  rate  sensor  data.  The  good  agreement  of  the  magnetometer-derived  pitch  rate  with  the  rate¬ 
sensor-measured  pitch  rates  whenever  those  measurements  exist  supports  the  accuracy  of  the 
magnetometer-derived  rate  estimates  at  all  other  times.  The  magnetometer-derived  rates  indicate 
initial  pitching  rates  of  approximately  4000  deg/s.  These  early  data  are  particularly  important 
because  the  high-drag  shape  of  the  model  causes  the  mach  3.5  launch  velocity  to  decay  to  subsonic 
speed  in  less  than  1  second,  and  the  high  mach  numbers  are  representative  of  re-entry  velocities. 
Later  in  the  flight,  as  the  model  begins  to  tumble,  pitch  rates  approaching  20000  deg/s  are  esti¬ 
mated  by  the  magnetometer  (figure  5b).  Although  these  data  are  not  of  interest  for  characterizing 
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CEV  aerodynamics,  they  demonstrate  that  the  magnetometer  method  does  not  suffer  from  dynamic 
range  limitations. 
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b)  Pitch  rates  during  flight 


Figure  5.  Pitch  rate  history  from  rate  sensors  and  magnetometers. 


With  6m ,  <f>m ,  q ,  and  r  in  hand,  equation  12  is  used  to  compute  y'/m  at  each  sampling  interval. 
The  magnetic  azimuth  is  then  given  by 


V m  (0  V m  (o) 


(14) 


Finally,  the  earth-fixed  Euler  angles  are  given  by 


re} 

(a  \ 

u  m 

V 

II 
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tq 

V  m 

K& s 

■jfrm  y 

(15) 


where  Tmf  is  the  DCM  relating  the  magnetic  and  earth-fixed  coordinate  systems.  This  metho¬ 
dology  has  been  successfully  implemented  in  an  on-board  digital  signal  processor  for  real-time 
guidance  of  experimental  projectiles,  as  reported  in  reference  1.  When  this  methodology  was 
executed  during  post-processing  of  the  flight  telemetry  data,  the  Apollo  model  heading  history  seen 
in  figure  6  was  computed  for  the  first  second  of  flight.  With  these  data,  aerodynamic  coefficients 
of  interest  were  estimated  for  the  test  vehicle,  and  the  flight  experiment  evaluation  was  successfully 
completed. 
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7.  Criteria  for  the  Use  of  Magnetometer-Based  Angular  Rate  Estimation 


The  effectiveness  of  this  methodology  is  clearly  dependent  on  the  accuracy  of  the  measurements 
of  M j  ,M  j  ,Mk ,  M  j ,  and  .  Thus,  a  calibrated  vector  magnetometer  is  required.  Calibration 

constants  can  be  determined  on  the  ground  and  pre-loaded  or  often  can  be  dynamically  detennined 
in  flight.  We  estimated  M j  and  M /(  by  differencing  successive  magnetometer  measurements. 

This  simplistic  method  requires  that  data  rates  be  sufficiently  high  to  accurately  estimate  the 
derivatives.  Sampling  rates  of  at  least  one  sample  per  degree  of  projectile  rotation  have  been 
found  to  be  adequate  for  a  number  of  simulated  projectiles.  Alternatively,  polynomial  fitting  to 
the  magnetometer  data  followed  by  analytic  differentiation  has  been  shown  to  produce  equally 
accurate  results  at  lower  sampling  rates.  Preferred  methods  should  be  determined  for  individual 
applications. 
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8.  Summary 


Free  flight  angular  dynamics  of  projectiles  have  been  successfully  measured  with  vector 
magnetometers  in  flight  experiments  during  intervals  when  angular  rate  sensors  have  failed  to 
provide  measurements.  This  result  argues  for  investigation  of  the  inclusion  of  magnetometers  as 
supplements  and/or  replacements  to  rate  sensors  in  low-cost  IMU/INS  systems. 
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